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Abstract. We use majorization and confidence intervals as a convenient tool to 
compare the uncertainty in photon number for different quantum light states. To 
this end majorization is formulated in terms of confidence intervals. As a suitable case 
study we apply this tool to sub- and super-Posissonian behavior and bunching and anti¬ 
bunching effects. We focus on the most significant classical and nonclassical states, 
such as Glauber coherent, thermal, photon number, and squeezed states. We show 
that majorization provides a more complete analysis that in some relevant situations 
contradicts the predictions of variance. 
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1. Introduction 

The most typical measure of uncertainty is variance, both in classical and quantum 
physics. However, variance is not the only measure of fluctuations nor the best behaved, 
as already pointed out for example in Refs. [T]. Some well-known disadvantages are, 
especially from a quantum perspective: (i) May fail when we move away from Gaussian- 
like statistics. In contrast to classical physics, non-Gaussian statistics is often the case 
in quantum physics, (ii) Puts too much emphasis on extremely unlikely outcomes. 
This may be a serious drawback in practical applications where the knowledge of the 
probability distribution is not complete, (iii) It cannot be used to derive Heisenberg-like 
uncertainty relations in finite-dimensional systems. 

This has motivated the introduction of several measures of quantum fluctuations 
alternative to variance EHS1E1EU7UH!. Most of them have the form of an entropy, such 
as the Renyi and Tsallis entropies, that include the Shannon entropy as a particular case. 
At the very heart of these alternative entropic measures lies the concept of majorization 


Sub-Poissonian and anti-bunching via majorization 


2 


since it is equivalent at once to the most popular entropic measures that satisfy the 
property of Schur-concavity. 

In this work we address the majorization of photon-number statistics in quantum 
optics. We put a particular emphasis on the comparison with Poissonian statistics. The 
reason for this is due to the fact that the comparison of variances is the most classic test 
disclosing nonclassical behavior in quantum optics and plays a relevant role in quantum 
metrology [9]. The most significant advances presented in this contribution are: 

(i) We show that majorization is equivalent to the unanimity of confidence intervals. 
Confidence intervals. They have been already used as a suitable measure of quantum 
uncertainty [Sj. This extends the equivalence of majorization to other measures of 
uncertainty beyond entropies. 

(ii) We show that majorization provides a much more stringent definition of sub- and 
super-Poissonian statistics, that may not coincide with the variance-based definition. 

(iii) We extend the analysis to cover bunching and anti-bunching effects. Using 
variance as a measure of comparison they are equivalent to super- and sub-Poissonian 
statistics, respectively. We introduce a suitable criterion for bunching and anti-bunching 
in terms of majorization as a more complete statistical assessment. We show that from 
this perspective the equivalence with super- and sub-Poissonian statistics no longer holds 
and may provide a new practical criterion to disclose nonclassical behavior. 

(iv) We apply this approach to the most frequent classical states, such as coherent 
and thermal states, as well as the most common nonclassical states, such as number 
states and squeezed states. All these states can be generated experimentally and have 
many practical applications. 

The deeper statistical characterization provided by majorization can have 
interesting physical consequences in quantum physics and quantum optics in particular 
since fluctuations and uncertainty are a cornerstone of the quantum theory. For the 
sake of illustration let us focus on possible physical consequences in the area of quantum 
metrology. 

One of the objectives of quantum metrology is to determine whether quantum 
fluctuations limit the precision in the detection of weak signals. Currently, most 
researchers believe that such limits exist under the form of Heisenberg limit pj. 
However, such belief is based on variance estimators, so alternative approaches to 
quantum uncertainty may eventually be useful to overcome their metrological effects, as 
already suggested in Ref. HU- The probe states with sub-Poissonian photon statistics 
are crucial to reach the ultimate quantum resolution. Therefore, a more complete 
understanding of sub-Poissonian statistics may shed light on these limits and their 
eventual breaking. 

Sub-Poissonian photon statistics is a paradigmatic nonclassical effect revealed by 
the comparison between variances. In this context, majorization is a comparison between 
complete distributions, so it may provide new tests of nonclassical behavior with a 
deep relation with previous nonclassical criteria introduced in terms of photon-number 
moments higher than variance dans]. Actually, majorization has been already used to 
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derive higher-order criteria of nonclasssical behavior in Ref. [ 13] . 

Deep down, we can regard quantum-metrology limits as the consequence of the 
uncertainty relations satisfied by complementarity observables, such as time and energy. 
Uncertainty relations are typically expressed in terms of variance, so it may be 
interesting to investigate the consequences that might be derived from uncertainty 
relations expressed in terms of generalized entropies, or even via majorization such 
as in Ref. |14j . For example, there is no position-momentum uncertainty relation if 
uncertainty is assessed using some particular entropic measures ra- Moreover, it has 
been shown that two observables can be complementary or not depending on the measure 
of fluctuations used [15] . and nontrivial uncertainty measures can solve long standing 
debates on the enforcement of complementarity 03. 

Finally, quantum metrology involves taking decisions with incomplete information. 
In such conditions it is preferable to work directly with probabilities rather than with 
their moments in order to obtain meaningful estimators of uncertainty. 

In Sec. 2 we recall the concept of majorization, relating it with confidence intervals 
as suitable uncertainty measures, as well as most common photon statistics. We compare 
different photon statistics in Sec. 3. In Sec. 4 we extend these ideas to bunching and 
anti-bunching phenomena. 

2. Majorization, uncertainty and photon statistics 

In this section we just recall the concept of majorization, its relation with uncertainty 
measures, and the photon statistics of the most usual and interesting classical and 
nonclassical states. 

2.1. Majorization 

The distribution p a majorizes the distribution p b , which will be expressed as pb -< p a , 
when the following relation between the ordered partial sums Sn holds, as illustrated 
in Fig. 1 

N N 

s n (pi) = E pij ^ 5Z itk = s n (pt), (i) 

j =0 k =0 

for N — 0,1... oo, where the superscript j. on /U denotes the same distribution p after 
a permutation of p n rearranging them in decreasing order 

Pi > pi > ■ ■ ■ > pj > ■ ■ ■ ■ (2) 

Majorization is a partial-ordering relation, so there are distributions that neither pb -< p a 
nor p a Pp b . 


2.2. Uncertainty, entropies and confidence intervals 

Roughly speaking, if p a majorizes p b we have that p b is more disordered or more mixed 
than p a . In other words, p f is more concentrated around its maximum. This can be made 
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Figure 1 . Relation between partial ordered sums Sjv and confidence intervals when 
the majorization pb -< p a holds. 


more precise by noting that majorization is equivalent to a relation between statistics 
in terms of doubly stochastic matrices [2j. This means that pb can be obtained from p a 
by the addition of classical noise. 

Thus majorization provides a very convenient measure of uncertainty. For example, 
the distribution with full certainty p\ = 1, p\^\ = 0 has Sn = 1 for all N, and majorizes 
all distributions. In the case of photo number, the states with full certainty are the 
photon-number states | n). 

Whenever two distributions can be compared via majorization, the result is 
respected by the Schur-concave entropies, such as Tsallis and Renyi entropies: 

R g(Pi) = yZ~ q ln ’ Tq ^ = T^~q (? P i “ X ) ’ ( 3 ) 

where q > 0 is the entropic index. The limiting case a —>• 1 is the Shannon entropy 
T\ — R\ — — J2i Pi log pi . Majorization implies that if pb -< p a then H q {j>b) > H q (p a ) for 
all q and H = T, R. 

As an alternative and more direct connection of majorization with uncertainty 
measures, let us consider the confidence intervals N(a). They are defined as the 
minimum number of indices such that the partial sum up to p^ N comprises a given 
fraction a of the probability 

N 

5»>cq N > N(a). (4) 

j=o 

When two distributions are comparable, pb -< p a is equivalent to saying that all 
confidence intervals of p\ are larger than or equal to those of p^\ 

Pb -< Pa <—» A T b (a) > N a (a) Va. (5) 

Otherwise, if the distributions are incomparable we will have Nb(a ) > N a (a ) and 
Nb(j3 ) < N a (f3 ) for some a, f3. 
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Let us recall the photon-number statistics of the most significant classical and 
nonclassical states beyond the trivial case of photon-number states. 


2.3.1. Coherent light The most celebrated light states regarded as classical in quantum 
optics are the Glauber coherent states. They have a Poissonian photon number 
distribution p c , n 


n 


Pc,n = —r exp (-n), 


A ~n = n. 


n\ 


( 6 ) 


This depends on a single parameter, the mean number of photons n, that equals the 
variance A 2 ?i. Unfortunately, there is no suitable expression for p\ n , nor for the partial 
sums SN(pj. n ). Nevertheless, in the case 1 the Poissonian distribution can be well 
approximated by a Gaussian, treating n as a continuous variable 

— (n — h) 2 /(2A 2 n) . (7) 


Pc,n - PG,n = 


An\Z2n 


exp 


In this case pQ n can be easily obtained because of the full symmetry and monotonic 
behavior around the maximum at n. The corresponding partial ordered sums are 


N 

Sn ( PG,n ) = X) PG,r, 


n =0 


2 1 V 

! V 2 An 

V^T JO 


dwexp w 2 ) . 


( 8 ) 


2.3.2. Thermal light Another basic and very common example of classical light are the 
thermal states, with photon-number statistics 

Pt^ = ~rr (~rr) » A 2 n = h(h+1), (9) 

n + 1 \n + 1/ 

that again depends on a single parameter. All thermal states are super-Poissonian 
A 2 n > n. The probability distribution is in decreasing order from the start p\ n = pt. n , 
and the partial ordered sums can be easily obtained analytically as 

Sn (ptn) = Y.pin = 1 - . ( 10 ) 

2.3.3. Squeezed light This is may be the best-known practical example of nonclassical 
light beyond photon-number states. Squeezed states are obtained from the vacuum 
|0) by a squeezing transformation S(r ) followed by a coherent displacement D(R) as 
D(R)S(r) |0). For simplicity we will consider squeezing parameters such that the photon- 
number distribution is [SJ 


Ps,n 


tanh”r 


exp 


— R 2 (1 — tanhr) 


x 


H„ 


2 n n\ cosh r 

tanh r- . 

y2 \ \J tanhr / 


( 11 ) 
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where R represents the coherent contribution while the squeezing parameter r, with 
sign, is the squeezed part that represents the squeezing of quadrature fluctuations. The 
mean value and variance of the photon number are 

n = R 2 + sinh 2 r, A 2 n = R 2 e 2r + - sinh 2 (2r). (12) 

They present sub- or super-Poissonian statistics depending on the R, r values. In any 
case there is no simple form neither for p\ n nor for the ordered partial sums Av(pin)- 

2-4■ Majorization applied to sub- and super-Posissonian statistics 

The idea of sub- and super-Posissonian statistics results from the comparison of the 
photon-number distribution p n of a given held state with a Poissonian statistics p c>n of 
the same mean number. This comparison has been always done in terms of variance. 
Majorization provides a much more stringent criterion as follows: 

p -< pc = under-Poissonian, p c -< p = over-Poissonian, (13) 

where in order to avoid confusions we have provided new different names. Over- 
Poissonian means that the curve S^(p^) is always above the curve S^ipj.) and vice 
versa for under-Poissonian. For Gaussian statistics over-Poissonian is equivalent to the 
sub-Poissonian and under-Poissonian equivalent to super-Poissonian. Otherwise, if the 
distributions are incomparable there can be no unambiguous claim about their relative 
photon-number uncertainty, and the result will depend on the particular entropy H q or 
confidence value a employed. 

3. Majorization relations 

Next we compare via majorization the photon-number distributions for the above 
classical and nonclassical states. When contrasting states from different families we 
will tend to compare states with the same mean number of photons. The reason to 
proceed in this way is that we are focusing on sub-Poissonian behavior, which implies 
comparison with a Poissonian distribution of the same mean. 

3.1. Coherent versus coherent 

It has been shown in Ref. [2] that any two Poisson distributions can be related through 
a doubly stochastic matrix in such a way that the distribution with smaller mean (this 
is smaller variance) majorizes the distributions with larger mean. Therefore, the order 
established by majorization coincides with the order established by variance and the 
distribution with smaller An always majorizes the distribution with larger An. As 
illustration in Fig. 2 we have represented the ordered partial sums Sn (pi,n) as functions 
of N for coherent states with different mean number of photons. This result is especially 
meaningful for small mean numbers. Otherwise, for n> 1 the Gaussian approximation 
holds and majorization readily follows from variance after Eq. (J8]) . 
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Figure 2. Ordered partial sums Sn (pi n ) as functions of N for coherent states with 
mean number of photons (from top to bottom) n = 1,5,10. Throughout lines are 
drawn as guides for the eye. 


3.2. Thermal versus thermal 


An equivalent majorization relation consistent with variance follows for thermal light 
directly from relation (10) 


Sn (pta) > Sn (p} lb ) <—» A ti0 n < A t , b n. 


(14) 


This case has been examined in more detail in Ref. |[I9]. In this case majorization may 
be pictured also as the result of added noise after Eq. (15) in Ref. [TS] for example. 


3.3. Thermal versus coherent 

Using variance as a measure of comparison, the thermal light should always have larger 
uncertainty than coherent light for the same mean number. However we have found that 
they can be incomparable for small mean photon numbers, such as for ii t = n c = 1.5, 
as shown in Fig. 3. The ordered partial sums cross at a critical value a — 0.80 for 
N(6t) — 3. Confidence intervals N(ot < a) say that thermal light has less number 
fluctuations than coherent light, i. e., behaves as over-Poissonian light, while the 
confidence intervals N(a > a) say the opposite and the thermal behaves as under- 
Poissonian. The transition occurs for a confidence interval that for Gaussian statistics 
would be between one and two standard deviations. On the other hand, for all cases 
examined with ii c = n t 1 the coherent states always majorize the thermal light, in 
accordance with variance. 

When we remove the condition of equal mean, we can get incomparable distributions 
even in the case of large mean numbers n c , fit 1. Moreover, in such a case we can 
find examples of full contradiction between majorization and variance, for example for 
h c = 100 and h t = 10. The variance of the thermal state is larger than the coherent, 
A= 110 versus A 2 c n = 100, while the thermal almost majorizes the coherent, as shown 
in Fig. 4. Strictly speaking, the ordered partial sums cross at a ~ 0.995 for N(a) — 57 
as shown in the inset of Fig. 4. However, the 99.5 % of the photon-number statistics, and 
thus almost all confidence intervals N(a), support that the thermal state presents lesser 
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Figure 3. Ordered partial sums 5% (pi >n ) and 5% [Pt, n J as functions of N for coherent 
(solid) and thermal (dashed) light with mean number of photons n t = n c = 1.5. 



Figure 4. Ordered partial sums Sn (pi, n ) and 5% ( Pt,nj as functions of N for coherent 
(solid) and thermal (dashed) light with mean number of photons n t = 10 and h c = 100. 
The inset shows an extension and detail of the plot for larger N revealing the crossing 
of the curves. 


uncertainty than coherent light, contradicting variance. This result clearly illustrates 
how variance overestimates outcomes with extremely small probability, such as 0.5% in 
our case. 


3-4■ Squeezed versus squeezed and coherent 


In the case of squeezed states we will consider always the same mean number of photons 
n for all states, but different number variance for some particular combinations of R 


and r Eq. (12). We have found incomparable statistics when contrasting sub-Poissonian 


states, when comparing super-Poissonian states, as well as when comparing them with 
Poissonian light. 

An example of incomparable states is illustrated in Fig. 5 where we have plotted 
the ordered partial sums for three different light states with h — 6 and variances 

Ag Ub n = 0.6 n (dashed), Ag Uper n = 14 n (dotted), and A 2 c n = n (solid). The super- 
Poissonian case is actually the squeezed vacuum R = 0. We can appreciate that both 
squeezed states are incomparable and the partial sums cross at d ~ 0.60 for N(a) — 2.6. 
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Figure 5. Ordered partial sums for variance-like super-Poissonian squeezed state with 
Aguper 71 = 14fi (dotted), sub-Poissonian squeezed state with Ag ub n = 0.6 n (dashed) 
and coherent light A 2 c n = n (solid) as functions of N with n = 6. The inset shows 
an extension and detail of the plot for larger N values revealing the crossing of the 
sub-Poissonian and coherent curves. 


Moreover the coherent and super-Poissonian states are incomparable as well and the 
partial sums cross at a = 0.85 for N(a) — 7. 

On the other hand, it seems in the main plot of Fig. 5 that the sub-Poissonian 
light (dashed line) majorizes the coherent light (solid line). Strictly speaking, this 
is not exactly the case if we look closer to large N values, displayed in the inset in 
Fig. 5, since the partial sums cross at a — 0.997 for N(a) — 14. This might be 
related to the characteristic photon-number oscillations of squeezed light. Nevertheless 
the fact that the states are incomparable has no practical consequences since it relies 
on extremely fine details of the photon-number distribution that will be hidden by 
experimental uncertainty. 

4. Bunching and anti-bunching 

A simple and direct practical observation of the sub- and super-Poissonian statistics is 
provided by the anti-bunching and bunching phenomena m- This is investigated by 
mixing the field state p with vacuum at a 50 % lossless beam splitter (see Fig. 6). At 
the output beams two photodetectors register the output photon numbers ri\ and n- 2 . 
In these conditions it can be easily seen that 



(15) 


Therefore, the sub- and super-Poissonian character is reflected unambiguously on the 
sign of AniAn^. super-Poissonian light produces positive correlations A 711 A 712 > 0, 
which is known as bunching, while sub-Poissonian light produces negative correlations 
AniAu 2 < 0, known as anti-bunching. 

We translate to this framework the ideas around majorization in order to provide 
statistical characterizations beyond An 1 An 2 . This turns out to be quite simple if we 
focus on the number sum n + — ri\ + 112 and number difference n_ = rq — U 2 since for 
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Figure 6. Arrangement for the observation of the bunching and anti-bunching effects. 


these variables we have 


AniAn2 = - (A 2 n+ — A 2 n_ 


(16) 


so that bunching is equivalent to A n + > A n_ while anti-bunching means An + < Arm. 
Because of number conservation and the vacuum input all registered photons ni, ri 2 
come from the state p and no photon is lost. Therefore, the number sum coincides with 
the number variable n + = n for all p. Then, after Eqs. (15) and (16) we have always 
A 2 n_ = n, so that the variable n + carries the number statistics of the state examined 
while the variable n_ provides a Poissonian reference. 

After Eq. (16) we can introduce a more complete criterion beyond variance via a 
majorization comparison between the distributions for the number sum p + and number 
difference p_. As we have just said, the p + distribution coincides with the number 
distribution of the input estate 

p+,n =Pn = (n\p\n), (17) 

while for we have 

1 


P-,m = Y, P- 


n=\m\ 



(18) 


and we note that when the number sum n is odd or even, the number difference m is 
also odd or even. Then we can introduce the following criterion 

p + -< p- = clustering, -< p + = anti-clustering, (19) 

where in order to avoid confusions with the variance-based bunching we will use the name 
clustering. Thus clustering means N + {a) > AL(a) for all a and is the majorization 
counterpart of bunching. Otherwise, if the distributions are incomparable there can 
be no unambiguous claim, and the result will depend on the particular entropy H q or 
confidence value a employed to assess the uncertainty. 

Note that the clustering criterion makes no reference whatsoever to the statistics of 
coherent states, since both p± in Eqs. ( flT| ) and (18) emerge exclusively from the number 
statistics of the state being examined. After Eq. (15) the variance-based bunching 
and anti-bunching are equivalent to super- and sub-Poissonian statistics, respectively. 
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N 


Figure 7. Ordered partial sums for the number-sum distribution p + (solid line) and 
the number difference distribution (dashed line) for a coherent state with n = 1. 


However this is not the case of the majorization-based clustering and anti-clustering 
with regard to under- and over-Poissonian statistics. Over-Poissonian means p c -< p 
while anti-clustering -< p. There seems to be no logical connection between these 
conditions since p_ in general is not Poissonian p c . 

Next we look for any relation between anti-clustering and nonclassicality that 
might parallel the well known role played by sub-Poissonian statistics as a nonclassical 
criterion. There is a simple connection if we focus on classical models where vacuum 
means a wave of zero amplitude. In such a case, after the 50% beam splitter in Fig. 6 the 
intensities recorded at both detectors are necessarily equal, regardless of the amplitude of 
the signal wave. This implies that the distribution has all probability concentrated 
in a single outcome, say m = 0, and therefore will majorize any other distribution. 
Thus, the anti-clustering criterion in Eq. (19) would be a nonclassical effect. 

As a first and readily example we may consider a coherent state. There is no simple 
expression for p_ >m in Eq. (18) so we resort to a numerical comparison. We have found 
that for large mean numbers, say n > 5, the partial sums of m and p\ j_ n are almost 
indistinguishable. On the other hand, for small mean numbers the situation is slightly 
different as shown in Fig. 7 for n = 1. The distributions turn out to be incomparable 
while all the partial sums of p£ m except the first one are below the sums of p\. n - The 
case h — 1 is where the difference of partial sums is maximal. Therefore, the classicality 
criterion p + -< p_ is not satisfied in this case, which otherwise is an example of classical 
light according to most of other criteria, with the exceptions in Ref. P2H This example 
also illustrates the lack of equivalence between clustering and Poissonian behavior. 

Furthermore, we have applied criteria (19) to the incomparable cases in Fig. 3 and 
Fig. 5 always obtaining majorization relations that agree with the comparison of their 
variances. We have obtained more interesting results examining the super-Poissonian 
squeezed state with n — 6 and A 2 n = 2 n. Despite variance, this state displays an 
effective over-Poissonian behavior (less uncertainty than a coherent state of the same 
n) for confidence intervals up to a ~ 0.9, as illustrated in Fig. 8. This agrees with 









Sub-Poissonian and anti-bunching via majorization 


12 



Figure 8. Ordered partial sums for the photon-number distribution of coherent 
squeezed state with super-Poissonian variance A 2 n = 2 n (dashed) and coherent light 
A 2 n = n (solid) as functions of N with n = 6. There is a clear effective over-Poissonian 
behavior (less uncertainty than a coherent state of the same n) for confidence intervals 
up to a ~ 0.9. This contradicts the comparison of their variances. 



Figure 9. Ordered partial sums for the number-sum distribution p + (solid line) and 
the number difference distribution p_ (dashed line) for the coherent squeezed state 
with super-Poissonian variance A 2 n = 2 n and n = 6. There is a clear effective anti¬ 
clustering effect (less uncertainty in the number sum than in the number difference) 
for confidence intervals up to a ~ 0.9. This contradicts the conclusions derived from 
variance. 


the ordered partial sums for p± in Fig. 9 that show that this state presents clear anti¬ 
clustering (less uncertainty in the number sum than in the number difference) for all 
confidence intervals with a values up to a = 0.9. 

In order to find further results contradicting the conclusions derived from variance 
let us consider the following mix of one-photon |1) and thermal p t states 

P = f|l)(l| + (! -Opt, l n )( n l’ ( 20 ) 

1 + n t ^ VI + ihJ 

where £ is a mixing parameter and n t is the mean number of the thermal state p t . The 
complete state p which has a mean number of photons of n — £ + (1 — £)hj. For example 
for n = 2 and £ = 0.9 the state p is clearly super-Poissonian A 2 n = llh. However 
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Figure 10. Ordered partial sums for the photon times thermal state (20) with n = 2 
and £ = 0.9 in dashed line and a coherent state of the same mean number in solid line. 



Figure 11. Ordered partial sums for the number-sum distribution p + (solid line) and 
the number difference distribution p_ for the photon times thermal state (201 with 
n = 2 and £ = 0.9. 


the whole statistics tell the full opposite. This is a consequence of the fact that the 
number distribution is strongly concentrated around n = 1, with p\ = 0.91, so for most 
practical purposes this works as a single-photon state, that is, over-Poissonian and anti¬ 
clustered. This is clearly reflected in Fig. 10 when comparing with a coherent state of 
the same mean number. The states are incomparable, but when referring to confidence 
intervals enclosing up to 90 % of the statistics a < 0.9 the state behaves effectively as 
over-Poissonian (less fluctuations than the coherent state) in spite of A 2 n = llh. This 
is again the effect on variance of large photon numbers with negligible probabilities. 

Regarding anti-clustering, the ordered partial sums of the distributions p± are 
represented in Fig. 11. They show that the states are incomparable. But again, for 
confidence intervals with a < 0.9 the state behaves as anti-clustered (less uncertainty in 
the number sum than in the number difference), in contradiction with the predictions 
of variance. In particular the probability that only one detector registers light —this is 
the probability of npn-i = 0 —is above 92%. 
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We have compared the photon statistics of the most relevant classical and nonclassical 
states from the perspective of majorization and confidence intervals. We have focused 
on the ideas of sub- and super-Poissonian statistics as well as the bunching and anti- 
bunching effects. 

As the most relevant results, we have found: (i) There are several interesting and 
common examples of incomparable states, (ii) Majorization may contradict variance, 
(iii) There is a majorization counterpart of bunching which is logically independent of 
Poissonian statistics. Let us briefly comment on these results. 

The information provided by majorization is valuable even when dealing with 
incomparable states. Previous works have found that different entropies lead to 
contradictory conclusions, where the minimum uncertainty states with one measure 
are the maximum uncertainty states of the other [22]. Majorization clearly explains this 
result showing that these contradictions arise because we are dealing with incomparable 
states [23]. 

Majorization can contradict variance because the former involves the complete 
statistics (that is, moments of all orders of the distribution), while variance only involves 
the first two moments. 

In any case, majorization will always provide a deeper knowledge of the situation 
than focusing exclusively on the variance. We do not propose to adopt majorization as 
a better way of comparing statics, but as another tool to understand quantum features. 
The use of one measure of the other will ultimately depend on the particular tasks 
addressed and the characteristic of the data available. 
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